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Verification of an acoustic transmission matrix analysis of 
sound propagation in a variable area auct without flow 

by J. H. Miles 

National Aeronautics and bpace Administration 
Lewis Research Center 
Cleveland, Ohio 4413b 
ABSIRACl 

a predicted standing wave pressure and phase angle profile for a hard 
wall rectangular duct with a region of converging-diverging area variation 
is compared to published experimental measurements in a study of sound 
propagation without flow. The factor ot 1/^ area variation used is of 
sufficient magnitude to produce large reflections. Ihe prediction is based 
on a transmission matrix approach developed tor the analysis ot sound propa- 
gation in a variable area duct with and without flow. 1 tie agreement between 
the measured anu preuicted results is shown to be excellent. 

NlIMt hCLAl Ukt 

matrix 

A area, m c 

c^ lsent^opic speed of sound, m/s 

t frequency, Hz 

H duct height equation parameter, in 

l unit matrix 

i i-n 1 ' 2 

k wave number, w/c , m * 

o o 

L 


L 


duct length, m 

length ot area variation, m 


m cl ln(A)/dx, m 

P transformation matrix 

p pressure. Pa 

k pressure reflection factor 

T transmission matrix 

t time, s 

u velocity, m/s 

x Cartesian coordinate, m 

V acoustic state vector 

y duct height, m 

l acoustic impedance, mks kayles 

ax region length, m 

x eigenvalue 

p q gas density, kg/m 

u angular frequency, radians/s 

Superscripts and Subscripts: 

( ) ' scaled quantity 

( reference state quantity 

( perturbed quantity 

INTRODUCTION 

An acoustic transmission matrix analysis of sound propagation in a 
variable area duct with and without flow is discussed in Ket . 1. Refer- 
ence 1 is concerned only with theoretical solutions of the sound propagation 
equations. This paper concentrates on the experimental verification of the 
analysis discussed in ket. 1. 
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The analysis is compared with an experimental standing wave pressure 
and phase angle profile obtained in a study ot sound propagation without 
flow in a long hard wall rectangular duct containing a region of converging- 
diverging area variation (Kef. if)- 

This experimental data provides a good test ot the theory since tne 
gradients of the area variation in the region of the converging-diverging 
area variation are large enough to produce a reflection much greater than 
the reflection produced at the duct exit. 

In the present paper, first a brief review of the theory is presented. 
Next, the experimental data is discussed. Then the analytical results are 
compared to the measurements and conclusions are drawn. 

THEORY 


The following assumptions are made, lhe duct has an area profile 
A(x). The acoustic equations are described by first order linearized, per- 
turbed equations of continuity momentum, energy and state tor isentropic 
sound propagation. 

State variable equation 

The perturbeo state variables are the pressure and velocity, lhese 
quantities are assumed to vary with time as exp (-iu>t). from Ret. 1, the 
Fourier transtorm ot the resulting equations is 
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Method ot solution 


Ihe method of solution used in this study is presented in kef. 1. for 
clarity, it is briefly discussed. The first step in obtaining a solution to 
tq. (1) is to divide the duct into a number ot regions or subsections. The 
region size is selected so that the area variation in each region can be 
approximated by an exponential area variation. Consequently, in a given 
region 


kith this assumption the value of n 1 n ( A) / ox is a constant in the region 
given by 

m =• d ln(A)/dx (3) 

Ihe b matrix in tq. (1) is now independent ot x in the region where it 
is to be evaluated. Thus , the solution to tq. (1) tor the jth region is 

f(Xj ♦ ax) = (e bAX ) j ^(Xj) (4) 


where x is the length ot the jth region ana exp (bax) is known as the 
matrix exponential. The value ot exp (bax) is oefineo by 


e bax 


1 + (bax) 



(bax) J 

— 


(S>) 


The matrix exp (bax) is the transmission matrix of the jth region (IK, 

J 

so that 



4 


Consequently, the transmission matrix, 1, for the variable area duct is 
founa from 


where 


Y(x) 


tt (bax) 


x=L 


Y(x) 


x=U 


T Y(x) 


x=U 


N 

c = Y2/ l AX )-j 
j=i J 


(?) 


18) 


Evaluation of matrix exponential 

While the series given in Eq. (5) defines the matrix exponential, it 
was not used for computer calculations since the convergence of the series 
may be slow. The matrix exponential was evaluated from 

e bAX = C" 1 P e bAX p _1 C (9) 

where 


D = p' 1 j/P (10) 

j* - Cbi' 1 (11) 


c = 


1 u 


0 p c 
o o. 
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P = 


' x i 


•**2 - 


•j/ 


n 


-si- 


ll 


( 13 ) 


and the a are the eigenvalues ot the si matrix given by the roots of 
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♦ bx ♦ c = 0 
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where 

b = -(j/ n ♦ j* n ) (lb) 

and 

c * ft j 2 ^ 22 ~ ^ 1 2 21 (1^) 

A more complete discussion of the evaluation of the matrix exponential can 
be found in Ref. 1. 

Boundary condition 

The standing wave profile is calculated as follows, first, the pres- 
sure at the duct exit (identified as location one) is assumed to be 1 Ra. 

A duct exit impedance is then used to calculate the velocity 

u l = 1 ll x ( 17 ) 


where the duct exit impedance is related to the duct exit pressure reflec- 
tion factor by 


7 / r 1 + k 
L l lo o c o = 


(IB) 


Next, the transfer matrices, T n , between each of the measurement points is 
calculated. The state vector at the last measurement point is 



T . T 
m-1 m-2 


... Wl 


1 

1 U 


1 


( 19 ) 


and the impedance at this point is 


' Pfi/ U m 


UU) 
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Note, that this is an initial value problem where the impedance and 
pressure are specified at the exit which differs from the boundary value 
approach used by two dimensional numerical solutions. 

The pressure at the point m-1 due to a unit pressure at the point m 
is obtained from 



( 21 ) 


This procedure is continued at each point. Finally, the pressure at the 
duct exit is 


pi 


- T~*T T - ^ 

- ij y-j ... 



( 22 ) 


The pressure magnitude o* the standing wave pressure profile is given by 



( 23 ) 


ana the phase angle of the standing wave phase profile is given by 


e 


tan 


jii^m\ 

k e p m J 


( 24 ) 


However, the pressure magnitude plotted is scaled with repect to the mea- 
sured pressure magnitude at the duct exit 



The phase angle is given with respect to the angle at x = 0. In addition 
since the experimental measurements of kef. 2 are maae assuming a time fac- 
tor of exp ( ♦ iu»t ) while the analysis is based on using exp (-iwt), the 
sign ot the angle calculated is reversed. Consequently, 

e « e(0) - e_ ( 2b ) 

m m 

THE EXPERIMENTAL DATA 

The general acoustic duct system described in Ref. 3 was modified into 
the simple no flow apparatus shown in Fig. 1(a) and used to measure the data 
discussed in Kef. 2. 

In the experiment discussed in Ref. 2 the duct temperature was 27 C so 
that the speed of sound was 348 m/s. The sound traveled in a 0.1 m by 
0.0381 m hard wall rectangular test section which was 1.42 m long. It was 
constructed so that the 8 top plates and 8 bottom plates were detachable. 

The variable area test piece was inserted in place ot the bottom piece at 
location 6. 

The variable area test piece is shown 'n Fig. 1(b). The surface pro- 
tile is prescribed by a fourth degree polynomial 

y/H = 16(x/L*)^ - 32(x/L*) 3 + 16(x/L*) 4 (27) 

where L* is U.U762 m and H is 0.0191 m which is just halt the duct height. 

The sound source used to obtain the data discussed in Ref. 2 was a 
speaker driven by a sine wave generator at lbbO Hz. The I860 Hz driving 
frequency assured that only plane waves would propagate in the straight por- 
tion of tne duct tar from the area variation since the first non plane mode 
begins propagating at a frequency ot 1740 Hz. An exponential horn was 
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attached to the upstream end of the test section. The horn was filled with 
an acoustic foam wedge to approximate a L * p Q c 0 impedance termination 
where the pressure reflection factor is K = 0.0 (Eq. (18)) which minimizes 
reflections ano simulates an infinitely long duct. 

Further information on the apparatus including a description of the 
microphone installation and the measurement procedure is given in Ref. 

DISCUSSION OF RESULTS 

The measured and theoretical normalized pressure magnitude profiles are 
compared in Fig. i and the measured and theoretical phase angle profiles are 
compared in Fig. 3. The theoretical curves calculated using as an exit 
boundary condition a duct exit pressure reflection factor of H = 0.0 which 
corresponds to an impedance termination of L * p Q c o are a near perfect 
match upstream of '.he area change, however, due to design and material 
limitations the pressure reflection factor is not exactly zero. Consequent- 
ly, the pressure magnitude is not constant downstream of the area change and 
the experimental ano theoretical curves are not in good agreement in tnis 
region. 

Impedance tube theory indicates that one unique value of the exit re- 
flection factor will generate the stanoing wave pattern shown in Fig. c. 

In order to determine this value, small changes in the exit boundary condi- 
tion were made so that the theoretical curves were in better agreement with 
the data downstream of the area change. The resulting boundary condition 
was R « -0.030 ♦ j 0.070 (Z * r p c p (0.94 •» j 0.0039)). This boundary 
condition as shown in Figs. ? and 3 yields a more perfect match not only for 
the downstream data points but at all data points. 
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The data discussed in Ref. 2 was used to verify the "steady" state 
f ,nite element theory of Ast ley and Eversman discussed in Ref. A ana the 
transient finite difference theory of White discussed in Ref. 5 for the 
limiting case of plane wave sound propagation without flow. Both of these 
numerical methods yield a two dimensional solution. It is shown in Ref. 2 
that both of these methods predict the proper trends of the data and are in 
excellent agreement with the measurements. In both cases, howe er, the cal- 
culations were made assuming that the impedance was l * p p c o . Con- 
seouently, the numerical calculations resembled the dashed curve in the 
downstream portion of duct (x greater than unity). 

CONCLUDING REMARKS 

The variable area transmission matrix theory (Ref. 1) producec good 
agreement with the axial data both in terms of the pressure magnitude and 
phe?o. Similar results can be obtained by numerical methods (Refs. 4 and 
5). However, these latter methocs are more complex than necessary if only 
plane wave propagation is of interest. Tho variaole area transmission 
matrix approach is simple to use hy comparison with two dimensional numeri- 
cal methods. Consequent ly, the acoustic transmission matrix analysis of 
sound propagation in a variable area duct is an ideal procecure to use for 
plane wave propaaation of sound in a hard wall duct with large area 
variat ion. 
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(a) No flow acoustic duct test section and instrumentation. 


(b) Variable area test piece with 5 microphone locations (microphone plugs not shown). 

Figure L - Apparatus. 





EXPERIMENTAL PRESSURE PROFILE 
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